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Abstract: We have explored Natural Supersymmetry (NSUSY) scenarios with low val¬ 
ues of the /r parameter which are characterised by higgsino-like Dark Matter (DM) and 
compressed spectra for the lightest MSSM particles, Xi; X2 xf- This scenario could 
be probed via monojet signatures, but as the signal-to-background ratio (S/B) is low we 
demonstrate that the 8 TeV LHC cannot obtain limits on the DM mass beyond those of 
LEP2. On the other hand, we have found, for the 13 TeV run of the LHC, that by op¬ 
timising kinematical cuts we can bring the S/B ratio up to the 5(3)% level which would 
allow the exclusion of the DM mass up to 200(250) GeV respectively, significantly extend¬ 
ing LEP2 limits. Moreover, we have found that LUX/XENONIT and LHC do play very 
complementary roles in exploring the parameter space of NSUSY, as the LHC has the 
capability to access regions where DM is quasi-degenerate with other higgsinos, which are 
challenging for direct detection experiments. 
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1 Introduction 

The naturalness of Supersymmetry (SUSY), which has been subjected to much discus¬ 
sion and thorough investigation for more than two decades [1-57], has become even more 
relevant now, as the Large Hadron Collider (LHC) collaborations ATLAS and CMS have 
started to probe SUSY in the TeV region. Indeed, the lack of evidence for superparticles at 
the CERN LHC, along with the rather high value of the Higgs boson mass in the context 
of the Minimal Supersymmetric Standard Model (MSSM), raises the question of whether 
the remaining allowed parameter space suffers from a high degree of fine-tuning, and if 
there is any parameter space of Natural SUSY (NSUSY) left. We discuss this problem in 
the framework of the well motivated MSSM. 

Based on standard measures of fine tuning [1, 2], the NSUSY parameter space was 
originally associated with light higgsinos, gluinos and stops, the SUSY partners of the SM 
Higgs, gluons and top quark. The present LHC limits on the masses of the latter two are 
approaching the TeV scale, under the assumption that the mass gap between these states 
and the lightest supersymmetric particle (LSP) , x?) is large enough (see for example [58- 
62]). Here Xi is the lightest of the four neutralinos, the mass eigenstates arising from 
the mixing of the fermionic component of the Higgs and gauge superfields, and 

(R^, W^), which are commonly known as the higgsinos, the bino and the wino. Moreover, 
these limits have a certain degree of model dependence and for stops, the experimental 
limits rely on certain decay channels being dominant (e.g. ii —)• tx)*, with t the top-quark) 
or on a substantial mass splitting between ii and x?) and can be significantly relaxed. For 
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example, in the scenario under consideration in this paper with higgsino like dark matter 
(DM), the stop branching ratios strongly depend on the left-right admixture of the lightest 
stop. Therefore model independent collider bounds on stops are weak or non-existent. 

It has however been shown that usual fine-tuning measures, defined as the sensitivity 
of the weak scale to fractional variations in the fundamental parameters of the theory, can 
be low even if the masses of the supersymmetric scalars are large. This happens in the so 
called “hyperbolic branch” (HB) [6] or “focus point” (FP) [14, 15, 63] regions of the minimal 
super gravity (mSUGRA) parameter space, where the value of the Higgs mass parameter, 
/X, can be low if the universal gaugino mass M 1/2 is not too large, as a consequence of 
the subtle interplay between the electroweak (EW) gauge couplings and the top-Yukawa 
coupling in the evolution of the squared Higgs mass parameters using the renormalization 
group equations (RGE). Moreover it was recently argued [64] that EW fine-tuning in SUSY 
scenarios can be grossly overestimated by neglecting additional terms, stemming from the 
ultra-violet (UV) completion of the model, that can lead to large cancellations favouring 
a low ^ parameter, but not necessarily a low stop mass up to a certain limit. Taking this 
point of view, we will take a low /x parameter to be the definition of NSUSY throughout 
our study. 

In the case /x <C Mi , M 2 (the EW gaugino mass parameters) one finds that the three 
lightest neutralino and chargino mass eigenstates, Xi^ X 2 xf, are quasi-degenerate 
and that these states are nearly pure higgsinos. In this scenario the DM relic density is 
typically below the WMAP [65] and PLANGK [66] measurements, because of the high rate 
of higgsino annihilation to standard model (SM) gauge and Higgs bosons and the higgsino 
coannihilation processes [67, 68]. One should note that this parameter space with the relic 
abundance below the experimental constraints is not however excluded, since the remaining 
relic abundance can be accounted for by other additional sources, e.g. axions. This NSUSY 
scenario, which is characterised by relatively light higgsinos in comparison to other SUSY 
particles, is not just motivated by its simplicity, but also by the lack of evidence for SUSY to 
date. We take advantage of the fact that NSUSY scenarios can be effectively described by 
a two dimensional parameter space, defined by the DM mass, i.e. the mass of x?) and AM, 
the mass difference between the DM candidate and the next to lightest supersymmetric 
(NLSP), typically xf, and our study explores the complementarity of the LHG and direct 
detection (DD) DM search experiments in covering this region. Such complementarity was 
the subject of recent studies, see e.g. [69] and references there in. 

It was already shown a decade ago that the HB/FP parameter space is challenging 
to probe at the LHG [70] even if the mass gap between Xi 2 xf is large enough to 
provide leptonic signatures. The most challenging case arises when the mass gap between 
these states is too small to produce any detectable leptons. The only way to probe such 
a scenario is via mono-object signatures, i.e. signatures involving a high transverse mo¬ 
mentum particle recoiling against missing transverse emerge (E™'®®), of which the monojet 
signature is particularly of relevance at the LHG, as initially suggested in [71] for generic 
compressed spectra. 

This technique has already been used in studies of quasi-degenerate higgsino spectra via 
monojetand monojet-|-E™®®-|-soft di-lepton signatures for the NSUSY parameter 
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space we consider [72-76]. However, we believe that these analyses are not entirely complete 
and/or have certain drawbacks. For example, in [72] the 95% confidence level (CL) reach 
for the 14 TeV run of the LHC was calculated assuming a signal (S) to background (B) 
ratio below the 2% level, which is probably not quite feasible when taking into account that 
the actual systematic error should be above 3-5% even for quite optimistic analyses [77-79]. 
In [74] the authors have performed their analysis at the parton level while, as shown in 
a preliminary analysis [80, 81], a fast detector simulation analysis leads to qualitatively 
different results and, therefore, is crucial. In [75] the conclusion about the observability 
of the quasi-degenerate higgsino NSUSY scenario from the monojet search was negative. 
However in this study the authors did not attempt to optimise the which turns out to 
be important as we will show in this paper. In [73] the prospects were more optimistic even 
after including systematics uncertainties. However, also there no optimisation of the cuts 
was performed which considerably enhances the accessible mass range as we will show. One 
should also note the Ref. [76] where authors suggested a new promising signature including 
a pair of soft leptons, and have demonstrated its potential power. However, in this study 
the important bb background was not considered as we discuss below, implying thus further 
background investigation for this signature which we do not consider at present. Finally one 
should mention Ref. [82], which studied similar to [76] monojet plus soft lepton signature 
suggesting visibly harder cuts to suppress bb background. The respective higgsino mass 
reach from this study is quite limited. On the other hand the suggested b-jet veto will 
not quite work for the signature with isolated soft muons, so, we believe that one should 
estimate bb background more precisely even for the case of harder cuts suggested in [82] . 

Motivated by the above-mentioned previous studies, here we aim to perform a com¬ 
prehensive and realistic analysis of the monojet potential to probe this NSUSY scenario. 
Our analysis is performed at the level of a fast detector simulation and the whole two- 
dimensional NSUSY parameter space mentioned above, rather than selected benchmarks 
as were attempted previously, is explored therefore completely covering the region of our 
interest. We consider prudent systematic errors and optimise the kinematic cuts to keep 
the S/B ratio at a reasonable level. We then discuss the LHC potential to cover the NSUSY 
parameter space at 8 TeV and produced projections for the 13 TeV run of the CERN ma¬ 
chine. By analysing also the exclusion potentiality for DM direct detection experiments, 
we aim to show that collider and DD experiments have a high degree of complementarity. 

The paper is organised as follows. In Section 2 we describe the parameter space and 
mass spectrum of NSUSY while in Section 3 we discuss the DM properties of this scenario. 
Section 4 is dedicated to the analysis of the collider phenomenology of the compressed 
higgsino scenario while in Section 5 we show the complementarity of collider and the DD 
experiment. We conclude in Section 6. 
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2 Parameter space and Spectrum of NSUSY 


In the bases {B^,W^, H^, H^) and the mass matrices of the neutralino and 

chargino sector of the MSSM are 


= 


/ Ml 0 —MzScoC/s MzSi^Sf^ ^ 

0 M 2 Mzc^^cp —MzCujSfi 

-Mzs^Ci3 Mzc^C 0 -fi 

\ MzSu)Sp —MzCuiSp —fi 0 y 


M-± = 

Xi 


/ M 2 \/2Mn/S/3 \ 




y 

( 2 . 1 ) 

where Mi and M 2 are the soft susy breaking mass parameter for B and W, fj, is the Higgsino 
mass parameter, Co; and Su, are cos and sin of the Weinberg angle, tan/I = V 2 /V 1 is the 
ratio of the vacuum expectation values of two Higgs doublets, Sjs, cp are sin/3 and cos/3, 
respectively and mz, mw are the masses of the SM gauge bosons and 

As a first step we consider scenarios where the Xi and have a high Higgsino com¬ 
ponent with all the SUSY partners of the SM fermions with masses in the multi TeV 
range and, to first obtain a qualitative understanding of the spectrum, we then expand the 
corresponding mass eigenvalues in the limit |/u| <C |Mi|, IM 2 I obtaining; 
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where we have defined S 2/3 = s\'n.{2(3)sgn{^) and a is the electromagnetic structure constant. 
In the case of xf we have also included the electromagnetic corrections as this shifts the 
mass by about 0.5 % which is indeed important in those cases where the mass splitting 
between xf/X 2 Xi is of the order of a few GeV. We would like to note, that in our 
numerical results we have included the complete one-loop corrections in the calculation of 
the masses as this changes the absolute mass scale by several per-cent. However, the mass 
splittings discussed below are hardly affected as the additional corrections changes those at 
most by 0(a/(47r) ln(m^/m^) which is below the per-mile level. For positive (negative) 
/i the mass eigenstate with negative (positive) CP-eigenvalue is the lightest one. The mass 
splittings are given by 


Amo = - m^o ~ 


Am± = m-± — m^o 

Xi Xi 


m| 


Ml M 2 


Amp ^ I Mmz) 

r% ' 


TT 


, rUy 
2 + hi^ 


(2.4) 

(2.5) 


where we have neglected corrections of the order 1/tan/3 and (^/Mi^ 2 )^. 

In order to analyse scenarios where the mass splitting between xf and Xi varies from 
the (quasi) degenerate regime up to the regime with larger mass splittings, we have chosen 
the following range of parameter space: 


fj, = (100, 300) GeV Mi = {fi, ij, + 600) GeV and (—//, —fi — 600) GeV tan/3 = (5 — 50), 

( 2 . 6 ) 
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fixing the value of M 2 = 2 TeV, which has the effect of decoupling X 4 and xf that will not 
be considered anymore in the following, along with the rest of SUSY spectrum, which is 
assumed to be decoupled. 

For I Ml I ~ /j, one can obtain a simple approximation for X123 masses [83] which we 
conhrm by numerical evaluation, indicating that in this parameter region x^ and X 3 are 
strongly mixed bino-higgsino states whereas X 2 is essentially a higgsino-like state with only 
a small bino component. 




0 50 100 150 200 250 300 

m^o [GeV] 




Figure 1. xf ~ Xi and X 2 ~ Xi niass splitting values as a function of m^o (upper row) and Mi 
(lower row) for the case Mi > 0 (left) and Mi < 0 (right). 

The mass splitting AM= m-± — m^o is shown in Fig. 1 as a function of and Mi 

Xi Xi 

for the case of positive and negative Mi while in Fig. 2 we show the contour lines for 
xf — Xi mass splitting. The relation between | Mi | — /r and the value of the mass splitting, 
which runs from quasi-degenerate scenario, AM~ 1 — 5 GeV for large Mi, to bigger values, 
AM~ 10 — 30 GeV for |Mi| ~ is clearly shown in this plots where the mass of the 
lightest neutralino, the DM candidate, is also presented. 

The nature of the neutralinos and charginos, as well as the small mass splitting between 
them, has a strong impact on their decay modes. In the case of pure higgsinos, the three 
body decays are dominated by virtual vector bosons, and due to the small mass differences 
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the decays into third generation fermions are suppressed. Note that in the scenario where 
Ml is close to |/r|, the off-shell lightest Higgs boson, hP, can also give sizeable contributions 
[84, 85]. Also, one should note that in the case |Mi| ~ /r, m^o — m^o ~ 2(m^o — 'm^o) and 
X 3 decays to the lightest chargino with the 50% probability while sharing about 25% decay 
to each of Xi and 


m^a and [GeV] 




Figure 2. Mass splitting between and (blue solid line) and mass of Xi (gray dashed line) in 
the fi-Mi plane. 


Three body decays in the limit of small mass separation are discussed in [ 86 ], where 
an effective theory study of the pseudo-Dirac DM scenario [87-89] such as the higgsino-like 
was performed. In this limit the decay width does not depend on the overall neutralino 
mass, but just on the mass difference 

^4 A 

r(x^x^^//'x?) = ^^^ (2.7) 

where A ~ is the mass of the leading mediator, while Am is either m^o — m^o or 

m-± — rn^o. As an example for off-shell Z exchange and decay into leptons, the coefficient 
C is 

+ ( 2 . 8 ) 

and a similar expression occurs for the case of the off-shell IT-decay. The proper decay 
length is very sensitive to the value of Am, and values below the GeV lead to displaced 
vertices, or collider-stable xf and X 2 - Indeed, for the decay X 2 //x? with an off shell 
Z exchange, the proper decay length is given by 

/ Am 

L = CT 0.025 cm f ^ ^ ^ j (2.9) 

which implies that for Am <0.1 GeV, X 2 would be collider stable. Similarly, for Am < 
1 GeV one could look for displaced vertices of order 100 fj,m. Note that the measured 


- 6 - 

















decay length would depend on the boost factor of the decaying neutralino to be taken into 
account as discussed in details in Ref. [86]. While collider stable X 2 contribute to the 
long lived or collider stable xf will provide a clear signature in the detector. A xf 
with a long enough lifetime can be detected in the tracking detectors by identifying decays 
that result in tracks with no associated hits in the outer region of the tracking system as 
recently analysed by ATLAS [90] and CM [91] collaborations. Both collaborations have 
obtained similar results, concluding on sensitivity for charginos with a lifetime between 0.1 
ns and 100 ns and covering chargino mass up to 500 GeV which significantly surpass the 
reach of the LEP experiments. For AM< 0.25 GeV the chargino could be a collider-stable 
charged particle [92], and bounds on such a situation arising from the 8 TeV run of the 
LHG can be estimated to be m^± > 300 GeV [93]. On the other hand, we have found 
that for AM> 0.4 GeV there is no limit on m~± > 100 GeV from the above LHG searches. 

Xl 

Therefore, our task is to analyse the potential of the monojet search to cover the NSUSY 
parameter space with AM> 0.4 GeV. 

3 Dark Matter direct and indirect detection in the NSUSY parameter 

space 

The results from Planck [66, 94] (see also WMAP [65]) have further decreased the error 
on the already very precise measurement of the dark matter relic density, = 

0.1184 ±0.0012. 

As we assume R-parity to hold, the LSP will be stable and will contribute to this relic 
density. In the scenarios under consideration, the LSP is the lightest neutralino, Xi, which 
is dominantly higgsino-like with a variable bino component. It is well known that for Xi of 
mass 0(100 GeV), a higgsino-like LSP provides DM relic density below the level observed 
by Planck. This happens because of the annihilation and co-annihilation rate of LSP and 
NLSP particles in the early Universe begin too large. On the other hand, for bino-like 
neutralinos the annihilation is suppressed, resulting in DM over-abundance. In the mass 
range we study, mixed bino-higgsino LSPs can therefore lead to the correct relic abundance. 
For /i < Ml the LSP is mainly higgsino-like and the value of ^ however 

this nonetheless at least solves the typical problem of the over-closure of the universe 
for neutralinos of this mass range in generic SUSY parameter space. In this case we then 
assume that the remaining relic abundance is accounted for by other means, for example, it 
could come from multi-TeV moduli field where the higgsino LSP is non-thermally produced 
(see e.g. Ref. [95]) or from mixed axion-higgsino DM (e.g. Ref. [30]). 

In order to assess the compatibility of the scenarios under our investigation with exist¬ 
ing experimental limits, we have evaluated the spin-independent annihilation cross 

section (ctsi) and the respective DD rates using micrOMEGAs 2.4.1 [96, 97]. In Fig. 3 we 
show the results for as a function of /x for tan /3 = 5, 50 and positive or negative val¬ 
ues of Ml. Note that the uncertainty on is not shown, the full one-loop corrections 

are not yet available, but we expect that these are not too large and will not qualitatively 
change our conclusions. From these plots we see, as expected, that in general lies 

below and decreases as Xi becomes increasingly higgsino-like. This is because 
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/Wi >0 /Wi <0 



50 100 150 200 250 300 50 100 150 200 250 300 


IJ [GeV] ij [GeV] 

Figure 3. The predicted value of the dark matter relic density is shown as a function of fj, 

for tan^ = 5,50 and positive (left) or negative (right) values of Mi as indicated. The relic density 
measured by the Planck satellite, is also shown for comparison, and the region excluded 

due to an overabundance of DM is indicated in grey. 


of the mass splitting between xf and x^, which becomes larger as Mi decreases. This 
suppresses the coannihilation channels which otherwise lead to an efficient reduction of the 
relic density. Below the WW threshold, the annihilation of the mixed gaugino-higgsino 
neutralinos therefore occurs via the Higgs and Z-bosons. The small bottom Yukawa cou¬ 
pling and the suppressed coupling to the Z implies that this mechanism is not efficient, 
apart from at the Z and h resonance. Therefore the spike in the relic density ^ ~ 100 GeV 
can be explained by the fact that this is just below the WW threshold. At lower values of 
H ~ 70, 90 GeV one moreover observes two dips corresponding to the Z-boson and Higgs 
funnels. This is most pronounced for positive values of Mi and lower values of tan/3 where 
the mass splitting is larger. 

In Fig. 4 we further show the spin-independent annihilation cross section for DD, again 
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Figure 4. The logarithm of the predicted value of the spin-independent annihilation cross section 
for DD (Tsi = Rn o'si/(10“® pb), rescaled by i?n/(10“® pb) where Rq = i® shown as a 

function of /r for tan /3 = 5, 50 and positive (left) or negative (right) values of Mi as indicated. The 
excluded limit from LUX (green), as well as the projected exclusion from XENONIT (cyan) are 
also shown for comparison, where the solid and dashed lines represent the exclusions for \Mi\ = fj, 
and I Ml I = /r-|-600 GeV respectively. 


for positive and negative Mi as in Fig. 3, where instead of crsi we plot the rescaled quantity 
Rn c^si (pb), where the scaling factor Rq = allows easy comparison with the 

most recent limits (also reported in these plots) from LUX [98], as well as the projected 
limits from XENONIT after 2 years live-time and 1 ton fiducial mass (see e.g. Ref. [99]), 
which in general assume the relic density to be the value measured by Planck. Fig. 4 
illustrates that the region with a low LSP masses and a higher mass splitting between 
xf and Xi (see Fig. 2), which is also easiest to see at colliders, is in fact excluded as 
the (Tsi for the DD experiments is too high. In the following we will further highlight 
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Figure 5. The mass splitting between xf/xi (blue solid line) and mass of Xi (black dashed line) 
is shown in the ^-Mi plane, along with the region excluded by LUX results (green). 

the interesting complementarity between the reach of the collider searches and the DD 
searches, particularly interesting for low DM masses. In Fig. 5 we also show the mass 
splitting between xf and x? in the (i-Mi plane of Fig. 2, along with the region excluded by 
LUX and the projected exlusion regions from XENONIT. This emphasises that the region 
where and Mi are very close each other is already excluded by LUX, which in turns puts 
an upper bound on the splitting between xf or x® and Xi- Note that for positive values 
of Ml, the splitting is larger, as is the mixing between bino and higgsino component of the 
DM. As larger mixing leads to larger couplings to Higgs bosons, DD is more sensitive to 
the case Mi > 0 . 

A related question is whether these scenarios could be excluded by indirect detection 
(ID) experiments, i.e. the detection of energetic e^, 7 , p or p, which may be created by 
the pair annihilation of weakly interacting massive particles (WIMPS). It turns out that 
the strongest bounds on neutralinos coming from such experiments are set by gamma ray 
telescopes: both the Fermi-LAT gamma-ray space telescope [100] as well as ground based 
telescopes. Fermi-LAT is sensitive to gamma rays particularly in the low mass range up to 
0(100 GeV). It is therefore particularly sensitive to lighter mixed bino-higgsino neutralinos, 
but the bounds are not competitive with those coming from DD. Here we are interested in 
light higgsinos with possibly some bino component, for which the relic density is in general 
much below the value measured at Planck/WMAP (see Fig. 3). In this case, all bounds 
from ID must be scaled by the square of the ratio of the predicted relic density to the 
experimental value. The rescaling appears as the pair annihilation cross section depends 
on the square of the local WIMP abundance. In [35, 101] it was shown that the Fermi-LAT 
limits derived for WIMP annihilations into WW (of which a large component of the total 
annihilation cross section should be comprised) are not yet sensitive to higgsino LSPs in 
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the 100-350 GeV range mainly due to the predicted under-abundance of neutralinos in 
these scenarios. Note that WIMPS having masses in the range 200 GeV to a few TeV 
will be probed by the future CTA array [102, 103], i.e. any higgsino or wino-like LSP for 
which the relic abundance is within an order of magnitude of Planck would be seen. In 
[104] a comprehensive scan of the pMSSM was carried out, and DM limits (both present 
and projected) coming from the LHG, DD and ID experiments were studied (note that 
the collider study was not dedicated to the region of interest of the present paper; our 
optimisation of the kinematical cuts goes beyond previous analyses). Ref. [104] found an 
impressive degree of complementarity: in the range a few hundred GeV to 1 TeV DD poses 
stronger bounds, and for higher masses ID is more sensitive. Therefore we conclude that 
for the mass range studied in this paper, the LHG and DD limits are the most relevant. 

4 LHC potential to probe NSUSY 

In the previous section we have discussed the current and future sensitivity of underground 
experiments to NSUSY scenario. In this section we explore the LHC potential to probe 
NSUSY and demonstrate that it plays a crucial complementary role. 

The LHC’s most sensitive searches for quasi-degenerate and Xi scenario are monojet 
signatures, i.e. the production of a pair of electroweakinos through the s-channel exchange 
of a SM EW gauge boson, Z, 7 or W^, accompanied by hard QCD initial state radiation 
(ISR) via the process 

PP -t XaXbj Xa,b = X?,2,3i xf ■ (4-1) 

The corresponding Feynman diagrams can be found in Fig. 6 . 



9 —3 ^ xt 9 kSiqsip -- 1 

Figure 6. Representative diagrams for pair neutralino-chargino production in association with 
quark/gluon leading to monojet signature. 

As we discuss below, the main problem for the signal search in this channel is the 
large background, the dominant contribution coming from the EW production processes 
Z-|-jets and IU-|-jets. We begin our analysis of the monojet signal from the recasting of the 
LHC at \/s = 8 TeV data and the respective experimental results, and then analyse the 
prospects of the 13 TeV LHC run with both standard and high luminosity (HL) options. 
Both ATLAS [78] and CMS [105] have performed studies of monojet signatures at the 
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LHC Runl, which have been interpreted in the context of an EFT approach. However 
this approach cannot be used for the NSUSY scenario because of the Z, W, 7 mediating 
interactions as indicated in Fig. 6 . Hence, in the case of NSUSY, 

DM with EW mediators =^Ceff = ^(^Fg) (xx) (4.2) 

where F is some Lorentz structure and x is the DM particle. Therefore it is necessary to 
recast the searches in terms of this NSUSY scenario and not to use limits from A-DMmass 
plane. 

4.1 Analysis Setup 

In this section we describe the different aspects of our simulation of both the signal and 
the most important backgrounds, implementing the important steps of hadronization and 
fast detector simulation. 

We performed a parton-level simulation using MadGraph vl.5.11 [106] with the MSSM 
model available on the FeynRules web page [107] implemented in UFO format [108]^, and 
cross-checked results against CalcHEP [109] with the MSSM model from the HEPMDB 
website^. 

At the level of matrix-element we have generated the production of a pair of elec- 
troweakinos via the s-channel exchange of a SM EW gauge bosons, accompanied by hard 
QCD initial state radiation. Parton level SM background simulations have been also cross 
checked between two packages. Our choice of PDF sets is CTEQ6L1 [110] and we used the 
MadGraph dynamical choice of renormalization scale which is equal to the geometric mean 
of Mass^ -|- P^ for the final state particles. Parton showering, hadronisation and decay 
of the unstable particles were simulated using PYTHIA v6.4 [111] while detector effects 
have been simulated with DelphesS [112] employing a suitable GMS card. Finally, the 
background processes yields, which include Z + j, W + j, tt, QGD and single top processes, 
have been taken from the experimental results for 8 TeV analysis and has been simulated 
using MadGraph-|-PTYHIA-|-DELPHES chain for the 13 TeV analysis. 

For the signal and background events at the parton level we have applied a cut on the 
jet transverse momentum of plp > 90 GeV at the generation level. This is a subtle point. 
One might argue that, since dealing with processes involving QGD radiation, it is necessary 
to apply a merging procedure between the hard jet generated via the parton-level matrix 
element and the soft jets generated by the showering algorithm, merging therefore the 0 -jet 
and 1-jet samples. However, since a hard jet is selected at the analysis level (e.g. plp > 110 
GeV for the 8 TeV LHG analysis and higher for the 13 TeV case, with a final selection 
requirement of a high somewhat correlated with the jet pT [80]) we found that this 

matching was unnecessary, and generated just the one jet sample. Moreover, avoiding 
matching in this case allowed us to make our analysis much more effective and obtain 
enough statistics in the high pj, region. To illustrate the validity of this procedure, in 
Fig. 7 we present the pT distribution of the leading jet for the case of X 1 X 2 production. 

^http://feynrules.irmp.ucl.ac.be/wiki/MSSM 
^http: //hepmdb. sot on. ac .uk/hepmdb: 0611.0028 
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Figure 7. Leading jet pt cross section distributions for the case of the 0+1 jet matched sample 
(blue) and 1 jet sample (black) for the X 1 X 2 production at the 8 TeV LHC. 


Shown in blue is the leading jet pp for the case of a 0+1 matched jet sample, while the 
case of the 1-jet events being unmatched is shown in black. One can indeed see that in 
case of a high pp cut on the leading jet at the analysis level, the matched and unmatched 
distributions of the leading jet are very similar, indicating that contribution from the 0-jet 
matched sample is negligible, and that the pp of the leading jet is dominated by the one-jet 
sample in the high pp region. 

4.2 LHC Runl: the reach of monojet searches 

We start our analysis with the exploration of the 8 TeV LHC potential to probe NSUSY 
and the recasting of the respective experimental results. For this purpose we have chosen 
CMS monojet analysis [105] which has been done for the data recorded at ^/s = 8 TeV with 
19.5 fb“^ integrated luminosity^. We have applied the following trigger selection followed 
by the cut-flow according to the CMS analysis: 

• Two triggers which require > 120 GeV or Fi™®® > 105 GeV and a jet with 

Pp > 80 GeV and within \r]\ < 2.6. 

• The analysis then requires that the jet with the highest transverse momenta has 
Pp > 110 GeV and \r]\ < 2.4. 

• Events with more than two jet with pp > 30 GeV and jr/j < 4.5 are discarded together 

with events where < 2.5, where ji and j 2 are the leading and sub-leading 

jets, to reduce QCD background. 

®The respective ATLAS analysis, which leads to similar results, can be found in [113]. 
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• The W production background was suppressed by applying a veto on events with 
one electron or muon satisfying pT > 10 GeV cut and events with one tau jet with 
Pt > 20 GeV and |r/| < 2.3. 

• Finally the analysis was performed in 7 regions with an increasing requirement of 

^miss ^ 250, 300, 350, 400, 450, 500 and 550 GeV. 

We have then derived the signal significance for each signal region from the number 
of signal events (S), after having imposed the cuts above, and the number of background 
events (B), using the following expression as in the GMS analysis: 

a = 2(^^STB - Vb) , (4.3) 


which is similar to the more common S/\/S + B for S' <C B, but is more robust with 
respect to downward fluctuations [114]. 



)U[GeV] 


Av[GeV] 


Figure 8. Contours for the 2(\/S + B — ^/B) (gray) and S/B (blue-dashed)in the plane (/r,Mi) 
for the signal regions 1 and 7 as defined in [105]. LUX and LEP exclusions are shown in green and 
red. 


In Fig. 8 we show contours of iso-significance a (solid gray) in the p-Mi plane for the 
lowest (left frame) and highest (right frame) requirements on The LUX exclusion is 

further shown by the green shaded area as well as the LEP2 limit on charginos by the red 
shaded area. In addition, the S/B ratio is shown by a blue-dashed line. By only inspecting 
significance contours, one could think that the Runl LHG data could slightly extend the 
LEP2 limit using the low E™®® signal region: the area on the left of the lines of a=2 would 
be ruled out at 95% confidence level (GL). However, one then observes that the S/B never 
goes above 2% in the region allowed by LEP2, while the actual systematic uncertainties of 
this analysis are of the order of 5-10% [78, 105]. This means that the low S/B ratio for the 
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8 TeV LHC is the main obstacle to going beyond the LEP2 limits. We can also see from 
Fig. 8 (right) that the higher cut increases the S/B ratio, eventually at the expense 

of the signal, meaning that the significance drops below 0=2 level in the parameter space 
allowed by LEP2. Therefore we conclude that Runl LHC does not have the potential to 
test the NSUSY seenario beyond the LEP2 limit due to systematie uncertainties. 


4.3 13 TeV LHC potential and complementarity to underground experiments 

As we have seen in the previous section, the LHC Runl is not sensitive to the NSUSY 
parameter space we consider due to the low statistical significance and the fact that the 
low S/B ratio remains below the systematic errors. 

In this section we study the LHC Run2 case and show that the higher collider lumi¬ 
nosity and energy allow us to choose kinematical cuts, bringing the S/B ratio to a desirable 
level while keeping the statistical significance at a high enough level in order to establish 
sensitivity to the NSUSY parameter space. Being a very challenging scenario, we make 
projections up to HL conhguration of the LHC machine for different assumptions regarding 
the E™*®® cuts and the control over systematic uncertainties. 

In order to ensure that the S/B ratio is under control, we compare the relative size and 
shape difference of the signal versus the dominant irreducible background Z+jet —>■ uv+jet 
{Zj). The relevant parton level distribution are shown as a function of the jet pT in Fig. 9. 
One can see that even for p = 93 GeV corresponding to m^o ~ 100 GeV, the background is 
about 3 orders of magnitude higher than the signal for a low plp > cut. An important feature 
of the signal versus background is that the shape of the background distribution is quite 
different from the signal: the background falls more rapidly with and the difference in 
the slope with respect to the signal is bigger for higher neutralino masses. The different 
in slope is mainly due to the mass difference between the neutralino, from the signal, and 
the neutrino, from the background. One should also notice that the difference between 
the shapes of the signal and background pj. distributions vanishes for very large values of 
plp » would expect. This provides us with our main strategy to optimise the 

LHC sensitivity to the NSUSY parameter space, that is to find the maximal value for E™'®® 
cut while maintaining the statistical significance at high enough level. On inspecting this 
figure one can already see that in order to achieve S/B ratio at about 5% level the value 
E™^®® cut should be around of 1 TeV. 

The signal process analysed is the same as that in Eq. (4.1), with the only differ¬ 
ence being that the initial requirement on the QCD ISR has been increased to 150 GeV, 
motivated by our preliminary study in Ref. [80]. 

Along with the signal processes, we have simulated the two main backgrounds 


PP Zj umj 
pp —>• Wj —)• I'^i'i -|- c.c. 


(4.4) 


with I = e, p,, T. 

Inspired by the 8 TeV CMS monojet search we have then applied the following cut-flow 
• We require a leading jet with px > 200 GeV and \ri\ < 2.4 
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Figure 9. Signal (dotted blue and dashed red) and Zj background (solid black) parton-level 
distributions for the 13 TeV LHC for the NSUSY scenario. Left: pip distributions for 100 fb“^ 
integrated luminosity. Right: normalised signal and Zj background distributions. 



Zj, Z —)• vv 

Wj, W 

p = 100 GeV 
Ml=700 GeV 

p = 200 GeV 
Ml=800 GeV 

Initial of events 

3.15-106 

1.25-10^ 

3.63-10® 

6.45-103 

pip > 200 GeV \r]^\ < 2.4 

1.05-106 

4.11-106 

1.73-10® 

3528 

Jet veto 

8.7-10^ 

3.13-10® 

1.33-10® 

2691 

A(/>(ji,j2) < 2.5 

7.2-10^ 

2.3 -10® 

1 .10-10® 

2320 

Veto ,pL^, 

7.2-10^ 

6.8 -10® 

1.08-10® 

2301 

> 200 GeV 

6.4-106 

4.3-10® 

9846 

2188 

> 600 GeV 

4353 

1002 

171 

93 

> 700 GeV 

1703 

250 

80 

47 

> 800 GeV 

694 

0 

37 

22 


Table 1. Outflow for the two main SM background and two choices of signal for the 13 TeV LHC 
with 100 fb“^ of integrated luminosity. The initial number of events corresponds to pip > 150 GeV 
cut. 


• We apply a veto on events with more than two jets with pp > 30 GeV and \rj\ <4.5 

• We require /S.(j{ji,j 2 ) <2.5 

• We apply a veto on electrons and muons with pp > 10 GeV. 

• We apply a veto on taus with pp > 20 GeV and \ri\ < 2.3. 

We have then defined signal regions with increasing cuts on an example cut flow is 

provided in Table 1. We have not generated tt, QCD and single top background, though 
we have applied to our simulated samples cuts that reduce these background to a negligible 
level with respect to Wj and Zj, see details in Refs. [78, 105]. 
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As observed for the case of 8 TeV LHC, a strong tension arises in attempting to 
simultaneous maximise the S/B ratio and the statistical significance. We demonstrate this 
in Fig. 10, where we plot S/B and a as a function on the final selection cut on the for 
both the case of x? ~ 100 GeV (left panel) and ~ 200 GeV (right panel). This figure clearly 
indicates that in order to achieve a high enough S/B ratio in keeping with the expected 
level of systematic uncertainties, a hard cut on ought to be applied, which at the 

same time pulls the significance down, due to the reduction in the number of signal events. 



Figure 10. S/B and a as a function on the final selection cut on the A™'®® for Xi ~ 100 GeV (left 
panel) and ~ 200 GeV (right panel). 

Given this tension it is therefore important to optimize the cut to provide a high 
enough S/B ratio and to keep o. >2 (5) in order to obtain an exclusion (discovery). The 
NSUSY parameter space can be conveniently described and presented in m^o-AM plane, 
where for AM we have chosen the mass difference between the lightest chargino and the 
DM particle. In this plane, for a given value of integrated luminosity, the optimal cut 
can be chosen by the point where S/B and a cross or are as close to each other as possible. 
This is related to the fact that the iso-significance contours are shifted to the left in the 
m^o-AM plane with the increase in the F^™®® cut due to the decrease of signal statistics, 
while iso-S/B contours are shifted to the right at the same time due to the increase of S/B 
ratio. Therefore the case when the respective iso-contours cross/are close to each other, 
would provide the maximal exclusion or discovery area in the m^o-AM plane. 

A1 

We illustrate this in Fig. 11 which presents S/B and significance isocontours in the 
m^o-AM plane for two different cuts on 850 and 900 GeV. One can see that indeed 

A1 

for F^™®® > 850 GeV, the exclusion area is below S/B = 3% (blue dashed) contour, while 
for FJ™®® > 900 GeV, the area below a = 2 (red solid) contour is excluded. Since for the 
first case the exclusion area is bigger, the F^™®® > 850 GeV is better choice for the optimal 
cut. We have found that a cut around 600 (850) GeV for 100 fb“^ (3000 fb“^) provides 
a ~ 2 and S/B 0.03 iso-contours optimally close to each other, which maximises the 
reach of the 13 TeV LHG for the NSUSY parameter space. The proximity of a ~ 2 and 
S/B ~ 0.05 iso-contours eventually requires higher F^™®® cut which is found to be around 
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LHC13 L=3000fb"^ 



/77^o[GeV] 

Figure 11. S/B (blue) and a (red) isocontours for two choices of if™'®® cut: 850 GeV (dashed) 
and 900 GeV (solid) in the m^o-AM plane. 

950 GeV, and as a result leads to a poorer 13 TeV LHC reach as we discuss below. While 
presenting results in the m^o-AM plane, we separate the cases Mi > 0 and Mi < 0, which 

Al 

differ due to different bino-higgsino components and mass splittings as discussed earlier, 

which eventually becomes less and less relevant as we approach the low AM region. 

In Fig. 12 and Fig. 13 we show results for Mi > /i and Mi < —/r (note that we have 

chosen /r > 0). The left panels contain the 2a exclusion LHC reach while the right panels 

contain the 5 ct discovery LHC potential. Both the cases of requiring 3% and 5% for the 

S/B ratio are shown, the latter just for the High Luminosity (HL)-LHC scenario. In the 

same plot we present also LUX sensitivity and the projected sensitivity of XENONIT to 

the in the m^o-AM plane. One can see that at Run2, with r-u 100 fb-i of integrated 
^1 

luminosity collected, the LHC will be able to exclude up to r-u 150 GeV Xi with AM below 
5 GeV, if the systematic error can be kept at the 3% level. The LHC will therefore surpass 
the LEP sensitivity for this scenario which reached the limit > 103 GeV for low AM 
values as shown in Fig. 12 and Fig. 13. By the end of the HL-LHC run, up to ~ 250 
(200) GeV LSP could be excluded for low mass splitting with a S/B> 3% (5%). It is quite 
remarkable that the LHC has the maximal sensitivity in the low AM region. This reach 
is nicely complemented by the LUX results and the projected exclusions for XENONIT, 
which cover the region with higher mass splitting. Such complementarity would allow LHC 
and DD experiments to completely exclude NSUSY scenario with DM mass up to ~ 250 
GeV. Discovery prospects for this scenario are shown in the right panel of Fig. 12 which 
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demonstrates that rsj 180 (110) GeV Xi can be discovered for S/B > 5% (3%) at the end 
of the HL-LHC run, while with 100 fb“^ the LHC will not have discovery sensitivity. 
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Figure 12. Exclusion (left) and discovery (right) contour lines for the 13 TeV LHC at the end of 
the LHC Run2 (light red region) and of the HL-LHC (light blue region) assuming S/B>3%. For 
the latter case also the case S/B>5% is shown. The region excluded by LUX and the projected 
exclusion by XENONIT are also shown, together with the LEP limit on the xf mass. Mi > is 
considered here. 

At the end of this section we would like to discuss a potentially promising new signa¬ 
ture involving a monojet plus soft di-leptons which was studied in terms of the NSUSY 
parameter space in Ref. [76]. In short, in this paper the signal from the second neutralino 
decaying leptonically into the lightest neutralino has been studied. It was suggested that 
one could trigger on soft leptons and use the upper cut on di-lepton invariant mass below 
10 GeV to suppress background and extract the signal from compressed Xi — X 2 produc¬ 
tion. After the suggested cuts, the background can be reduced down to about the 6 fb 
level bringing it below the signal, which would allow one to claim a discovery for chosen 
benchmarks with only 100 fb“^ of integrated luminosity. One should note however, that in 
spite of the comprehensive set of backgrounds, jbb was not considered in this paper. After 
application of the cuts from [76] we have found that jbb background is in fact dominant 
for the monojet plus soft di-lepton signal and is about two orders of magnitude above the 
backgrounds taken into account in the aforementioned study. An estimation of jbb back¬ 
ground is not trivial, since the cross section of the jbb process is very high, about 100 nb 
for soft initial cuts on jet pT, while the efficiency of the selection cuts is very low, about 
10“^. Therefore, a reasonable estimate of this background requires either the simulation 
of at least 10® events or the extraction of this background from future experimental data. 
Our preliminary results indicate that one should develop a dedicated strategy to suppress 
this jbb background in order to make the monojet plus soft di-lepton signal a viable tool 
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Figure 13. Exclusion (left) and discovery (right) contour lines for the 13 TeV LHC at the end of 
the LHC Run2 (light red region) and of the HL-LHC (light blue region) assuming S/B>3%. For 
the latter case also the case S/B>5% is shown. The region excluded by LUX and the projected 
exclusion by XENONIT are also shown, together with the LEP limit on the xf mass. Mi < —fj. is 
considered here. 

for the exploration of the NSUSY parameter space. 

5 Conclusions 

In this paper we have explored the complementary potential of the Large Hadron Col¬ 
lider and underground experiments to probe Dark Matter (DM) in the Natural Super- 
symmetry (NSUSY) scenario. This study, which combines searches from different kinds 
of experiments, has to be done in the context of a specific model, as (model-independent) 
Effective Theory (EFT) approaches are very limited in scope, see e.g. the discussion in 
Refs. [115, 116]. In particular the EFT approach is not applicable for well motivated 
NSUSY scenario, which we study here, where DM has direct couplings to Standard Model 
electroweak (EW) gauge bosons and the Higgs. 

Current limits on simple SUSY scenarios are at the TeV range, in clear tension with 
naturalness arguments and hence with the motivation for introducing SUSY in the first 
place. A possible explanation for this situation is that the manifestation of SUSY is not 
as simple as one expects, but there is more complexity in the structure of SUSY at high- 
energies. Notwithstanding, one would still expect that the particles more directly related 
to the tuning of the EW scale remain light in the spectrum. This leads to a generic 
expectation that DM in NSUSY should have a sizeable Higgsino component. 

While being theoretically attractive this scenario also represents a clear example of 
how colliders and underground experiments can complement each other. Indeed, while 
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underground experiments have a larger mass sensitivity than colliders, being able to probe 
the multi-TeV region, colliders can cover parameter space hidden from DM direct detection 
(DD) experiments. Specifically, the increase of the DM higgsino component makes NSUSY 
parameter space increasingly difficult to probe in DD experiments. In this region higgsino- 
like DM is quasi-degenerate with two other particles, the second neutralino and the lightest 
chargino, and the increase of higgsino component is correlated with the decrease of this 
mass splitting, AM. At the same time, the LHC sensitivity increases with the increase 
of the higgsino component of DM and reaches its maximum for very low values of AM. 
We have conveniently described the NSUSY parameter space in the fi-Mi region and have 
translated it into sensitivity of the LHC and DM DD experiments in the m^o-AM plane, 
where for AM we have chosen the mass difference between the lightest chargino and DM. 
We have studied the current and the future reach of underground experiments in this 
region, as well as given predictions for the relic abundance in NSUSY. 

We would like to point out that we have presented combined LHC and DM DD exper¬ 
iments results for the whole m^o-AM NSUSY space, rather than for chosen benchmarks. 
Moreover, we have optimised the final selection cut to keep the S/B ratio at a high 

enough level, to deal with the systematic errors on the background, which aer one of the 
main problem for the exploration of the NSUSY parameter space. As a result, we have 
found that the 8 TeV LHC unfortunately does not allow to test the NSUSY parameter 
space beyond the LEP2 limits, while the 13 TeV LHC at 3 ab“^ has the potential to 
significantly surpass the LEP2 reach and to cover DM masses up to about 250 GeV for 
AM<5 GeV, which can neither be covered by LUX nor by XENONIT DM DD experi¬ 
ments. At the same time the XENONIT experiment will be able to complementarily cover 
the AM>5 GeV parameter space up to large values of the DM mass, well beyond the LHG 
reach via mono jet analyses. 
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